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Abstract. We present the phase structure of the chiral 577(2) x SU(2) scalar- 
fermion model on the lattice using the Zaragoza proposal for chiral fermions. The nu- 
merical result agrees with an analytic study based on the use of weak and strong yukawa 
coupling expansions combined with the mean-field approach. The phase diagram con- 
sists of four phases: paramagnetic (PM), ferromagnetic (FM), antiferromagnetic (AFM) 
and ferrimagnetic (FI). The transition lines separating these four phases intersect at one 
quadruple point. 
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1. Introduction 



The formulation of a Chiral Gauge Theory (CGT) on the lattice suffers from the 
well-known doubling problem [1]. Several ways of dealing with this problem have been 
reviewed in reference [2]. Among them, the Zaragoza proposal [3,4] is a viable candidate 
for describing chiral fermions on the lattice. This proposal belongs to a category of 
models wherein, as in the Rome approach [5] , gauge symmetry is explicity broken by the 
regularization and therefore additional counterterms must be included in order to restore 
the gauge invariance in the continuum limit. The distinguishing feature of this approach 
is that the global chiral symmetry is preserved, and the doubler fermions are free and 
massless. These properties simplify the calculation of some important parameters of the 
Standard Model (the S 7 U and Ap parameters) [6] . 

In the case of Chiral Yukawa models (CYM), this regularization method preserves 
the invariance under discrete rotations and translations, but also all the important sym- 
metries present in the continuum: hermiticity and global chiral symmetry. Hence it is 
particularly well suited to study those models [4] . 

In this paper, we investigate the phase structure of a SU(2) L x SU(2) R CYM which 
is essentially the fermion-scalar sector of the electroweak theory. We decided to freeze the 
radial mode of the scalar field, which corresponds to the choice of an infinite bare quartic 
coupling. Knowledge from the pure <& 4 theory suggests that such a model belongs to the 
same universality class as the models with finite quartic coupling. For a comparison of 
the phase structure of different models see references [2,7-10]. 



2. Model, symmetries and limiting case 

Our action S(^ 7 $) is given by 
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(1) 



where 

M*) = ~\ E Tr (^A** + ( 2 ) 
is the kinetic term for the scalar fields, <&, which are 2x2 SU(2) matrices. 

S f(^) = \J2 (^xT/i^+A - Vx+itf^x) (3) 

is the fermionic kinetic term where J2y stands for the summation over rif doublets of 
Dirac fermions, and finally, 

s Y (v, *) = y £ + (4) 

is the Yukawa interaction. The L and i? indices refer to the left and right components of 
the fermion fields. In this interaction term, the way to implement the decoupling of the 
doubler fermions is based on the use of the quasi-local field component given by 

x' £ hc(x) 

Here he (x) is the elementary hypercube that starts from the site x in the positive direc- 
tion. To understand why we have chosen to describe the fermion interactions in Eq. 
(4) let us go to momentum space. There is given by 

& 1 Hq) = F(qMq), F(q) = ]Jf( q ,), f(q„) = cos(|), qE(-7T^} d . (6) 

Note that is the original field \I> modulated by a form factor F which is responsible 
for the decoupling of the doublers at the tree level in the continuum limit. This is achieved 
by forcing the form factor to vanish for momenta corresponding to the doublers (note 
that /(0) = 1 and f(ir) = 0). Other form factors are possible but eq. (6) corresponds to 
the most local possible choice for the smearing in position space, given by eq. (5). We 
have proved perturbatively [3,4] that in the continuum limit only one fermion (for each 
flavour) is coupled to the scalar. The doublers stay free and massless. 

The action (1) has a global chiral symmetry and is also invariant under a shift 
transformation of the fermion fields. This shift- symmetry can be used to argue, if a 
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continuum limit exists, the non-perturbative decoupling of the doublers (see Ref. 4 for 
more details). 

A few words are needed about the symmetries of the phase diagram because, in 
this model, they are not exactly the same as in other regularizations and have relevant 
consequences. 

Some of the usual symmetries are present. The action is invariant under the trans- 
formation & x — > — y — >• — y, so we can restrict ourselves to y > with no loss 
of generality. For y = the action is invariant under k — >■ —k, <& x — > e x & Xj e x = 
^_^x 1 +x 2 +x 3 +x 4 ^ ^ e scg j ar an( j f erm i on i c fields are decoupled and the model reduces to 

a <& 4 model with the radial mode frozen: two critical points are present at k = ±k c . 

But the symmetry k — >• — k, y — >■ —iy, usual in other models, is absent here [4]. 
As a consequence the phase transition line PM-AFM is not determined by the phase 
transition line PM-FM, as often happens [7]. The simple mean field computation that 
we shall present in sect. 3 is able to detect the crossing of these lines and the presence 
of a ferrimagnetic phase. 

Concerning the limiting case y — >• oo, we have a peculiarity similar to that in the 
Chiral Yukawa model with hypercubic coupling [9]. A rescaling of the fermionic field 
\I/ — >• (l/y 7 ^) suppresses, as always, the fermion kinetic term but, in our case, the 
fermionic fields on different sites are still coupled by the Yukawa interaction Sy (because 
of the smearing in ty^). The fermion can propagate and does not decouple from the 
scalar field when y goes to oo. In this limit the model is not equivalent to a pure <& 4 
model. 



3. Yukawa coupling expansions and the mean field technique 

To evaluate the free energy, we have used conventional mean-field methods for the 
scalar field [11] together with weak and strong coupling expansions for the Yukawa inter- 



action. At small y, the Yukawa term: 



exp 



(7) 



is straightforwardly expanded in powers of ?/ up to the four fermion interaction (i.e. up 
to y 2 ). 

On the other hand, for large values of y, we first introduce auxiliary fermionic fields 
rj x and rj x using the following identity: 

y 2d/2 J dr} x dr\ x exp j l - (rj x & x P R r} x + rj x $ x P L r} x ) + ^ { x \ x + rj x *W J , (8) 

and proceed with an expansion of the exponential term in powers oil/y. Here again, we 
have kept all the terms up to the four fermion interaction (i.e. up to 1/y 2 )- 

The next step is to make the scalar field dependence of Sb(&) linear by using the 
auxiliary fields V x , V^, A x , A x (N is the number of sites): 



exp 



I x,n ) 



( — \ / [dVdVUAdrf] exp - Tr (t£ V x+A + V^V Q 



+ iE^[4(^-^)+4($i-vi)] 



(9) 



In both cases (for small or large y ), the functional integration over the scalar field 4> 
can be performed and the fermion terms exponentiated. Before doing the integration over 
the fermion fields (in order to convert the Yukawa model into a purely bosonic model), 
we need to decouple the composite fields of the type ^ x ^ x (which appear in the four 
fermion interaction) by using identities like [12]: 

expJijE*^ 1 ') V 1 } = / ^expJ-iA^ + M-^X:*! 1 '*^} (10) 



As a result of the bilinear structure of the action in the fermionic fields, it is now 
possible to carry out the fermionic integrations in the partition function. The remaining 
integrals are approximated by the mean field saddle point method: we look for a transla- 
tionally invariant saddle point with a constant and a staggered piece, we then substitute 
the auxiliary fields V x and with v + v st e x , A x and A\. with —i(a + a st e x ) and A^ with 
X + X st e x . This approach yields the free energy per unit volume, 

1 . _ 1 



N ° 2 



T = log Z = ^(X 2 + X 2 t )-2k(v 2 -v 2 t ) + av + a st v st 

^ [u(a + a st ) + u(a - a st )} - ^2 d / 2 /, 



(11) 



with k = kd 7 u(/3) = log ^Ii(f3), where I\ is the modified Bessel function of order 1. The 
function / comes from the fermionic determinant and has different expressions for strong 
and weak yukawa coupling, I = Is and / = Iw respectively with: 



Iw= J ^logj[ S 2 (p)+^^^ 



(12) 
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Is = J ^ d log^[zSs 2 (p) + y 2 F 2 (p)F 2 M] 2 + y 2 ^j s 2 (p)(F 2 (p) - F 2 (p)Y 

(13) 

where F 7T (p) = F(p x + tt, . . . ,p d + ir), and 

s 2 (p) = ^2sin 2 p x , 



z = u(a + a st ) - (A + X st )\/u(a + a st ), 



z = u(a - a st ) - (A - X st )\/u{a - a st ). (14) 

The saddle point equations are obtained by requiring that the free energy is extremal 
for v, a, A, v st , a st and X st . The free energy in the presence of source terms can be evalu- 
ated following the same steps, and the mean field predictions for the order parameters are 
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obtained by taking derivatives with respect to these sources. For instance, for a lattice 
with N sites, with the definitions 

<*> = (^£**)> <*-> = (^ ( 15 ) 

we find (<&) = f , ($ s t) = f s t if we align the symmetry breaking direction along the unity 
matrix. Looking numerically for the saddle point values and the free energy for several 
choices for the parameters k and y, we have found four phases, (see fig. 1): 
Paramagnetic (PM): Here, (<&) = (<& st ) = 0. 
Ferromagnetic (FM): In this phase (<£>) 7^ and (<& s t) = 0. 

Antiferromagnetic (AFM): This phase is characterized by (<3>) = and (<fr s t) / 0. 
Ferrimagnetic (FI): Both parameters different from zero. 

The expansion of the term e~ Sv is essentially an expansion in y ^f^^f^ or in 
1/y rj x rj Xl and as a consequence we expect that our results are not too bad in regions 



w 



here 



y ( ty^^x^ ) I < 1 at small y or where 



(1/y) (rg^ 



< 1 at large y. In the 



figure we have only shown the regions of y where these conditions are fulfilled. 

One last comment about the mass of the physical fermions in the FM phase, mp 
(in lattice units). For small y, it can be shown, from eqs. (11,12), that m F = y 2 z 2 . 
For large y 7 eqs. (11,13) gives mp = y 2 / z 2 . Therefore, as usual [7], we can distinguish a 
weak and a strong FM regions. 



4. Phase transition lines 

A sampling of the values of the order parameters for several values of y suggests 
that the change in the order parameters occurs continuously when crossing the phase 
transitions. Close to the phase transitions, this allows a linearization of the saddle point 
equations in the mean field variables which vanish at the phase transition (those are v, 
a, A for the FM-PM and AFM-FI phase transitions and v 8t , a st , X st for the AFM- 
PM and FM-FI phase transitions). In the phase diagram plotted in the figure, the 
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FM-FI and AFM-FI transition lines have been obtained with a numerical solution of 
the linearized equations. For the FM-PM and the PM-AFM transitions an analytical 
solution is possible; the transitions are of second order and the critical lines are: 



where the upper signs stand for the FM-PM transition and the lower signs for the PM- 
AFM. For d = 4, 



Because of the very small value of /_ the PM-AFM transition line is almost straight. 
One would expect this behaviour in any theory in which the doubler fermions are decoupled. 
In fact, in the mean field approximation, without any form factor, in the AFM phase, 
the Yukawa term couples the scalar mean field, v s t, to the fermion sector through a 
coupling of the type ^ k'Vst^ kO~(k — k' — tt). Because of the factor S(k — k' — tt), the only 
coupling between the scalar and the fermion is through a vertex involving simultaneously 
the physical and the doubler fermions. Therefore, if somehow the doubler fermions are 
effectively decoupled from the scalar sector, simultaneously the scalar sector is decoupled 
from the physical fermion, and the PM-AFM phase transition line will be straight. 

In our model it is easy to understand how this phenomenon occurs. In I + , F 4 (9) 
admits the contribution of the physical fermions and kills the contribution of the doublers, 
whereas in F 2 (6)F 2 (0) kills the contribution of the doublers but also that of the 
physical fermions. The PM-AFM transition is nearly independent of y. We can compare 
with the case of the naive fermions where F 4 (9) = F 2 (9)F 2 (9) = 1, so i+ = I- and 
we recover the k — > —k, y — > —iy symmetry; the transition lines PM-FM and PM-AFM 
are now parallel, the FI phase and the quadruple point are not seen in the mean field 
appr oximat ion . 

Note also that the average slope of the PM-FM transition line in our case is smaller 
than when the doublers are coupled. Consequently the value of the bare Yukawa coupling 



(16) 




(17) 




(18) 
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at the quadruple point is large. Moreover, perturbative renormalisation group arguments 
are in favor of an increase of the renormalized coupling constant when the number of 
coupled fermions decreases. This could be relevant for the upper bounds of the physical 
masses. 



5. Monte-Carlo results 



We have performed an exploratory numerical simulation of this model on lattices of 
size 4 4 and 8 4 . We have used the Hybrid Monte Carlo algorithm [13] with two flavour 
doublets (an even number is required for the algorithm). The fermion matrices was 
inverted with the Conjuguate Gradient algorithm. The length of the molecular dynamics 
trajectories was selected with a random distribution with a mean value of 10 steps. For 
chosen values of y, we scan over A; in a region where a phase transition is predicted by the 
mean field results given in the previous section. For each simulation in the (y, k) plane, 
the number of configurations varies between 3000 and 14000 after having discarded about 
1000 trajectories for thermalisation. We have measured the order parameters and the 
associated susceptibilities. In the large y region, we have only checked that we have 
a ferrimagnetic phase when k is sufficiently negative without determining precisely the 
transition point between the FM and the FI phases. In the small y region, the measured 
transition points are plotted on the figure and compare well with the results from the 
mean field predictions given by the dashed lines. 



6. Conclusions 

We have determined the phase structure of an SU(2) x 577(2) fermion-scalar model 
with the globally chiral invariant model proposed in [3,4] to decouple the doublers. The 
numerical and mean field results are in good agreement. We have found four phases (PM, 
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FM, AFM, FI) with a quadruple point. The mean field calculation is able to predict the 
FI phase and the quadruple point. The PM-AFM transition line is nearly independent 
of y, this property is intimately related with the decoupling of the doublers. 

The global chiral symmetry is preserved by this model and the decoupling mechanism 
allows the PM-FM transition line to reach rather large bare yukawa couplings. Thus, 
the questions of a bound on the mass of a heavy fermion induced by a strong Yukawa 
interaction and its influence on the Higgs mass should be fruitfully investigated along 
this way. These issues are currently under investigation. 
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Figure Caption 

Fig.l Phase diagram for the SU(2) L x SU{2) R Chiral Yukawa model, for rif = 2 doublets 
of Zaragoza fermions. The dashed lines are the mean field results for the transition 
lines. The open symbols are the transition points determined by the Monte Carlo 
simulation. The square correspond to the FM-PM or AFM-FI transitions, the circles 
to the AF-FM or FI-FM ones. Because of CPU limitation, we have not succeeded 
in a precise determination of the position of the quadruple point. 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-lat/9402001v2 



